A Levinson theorem for scattering from a Bose-Einstein condensate 
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A relation between the number of bound collective excitations of an atomic Bose-Einstein con- 
densate and the phase shift of elastically scattered atoms is derived. Within the Bogoliubov model 
of a weakly interacting Bose gas this relation is exact and generalises Levinson's theorem. Specific 
features of the Bogoliubov model such as complex-energy and continuum bound states are discussed 
and a numerical example is given. 
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Exact results in physics are rare. One such result is 
Levinson's theorem yj, relating the number of bound 
states of a given potential to the accumulated scattering 
phase shift at threshold. This theorem was first formu- 
lated for the single-particle Schrodinger equation and in 
general does not hold for many-body systems. Instead, a 
generalisation for multi-channel scattering exists 0], re- 
lating the cumulative sum of the phase shifts of all open 
channels to the number of bound states. Little can be 
said about phase shifts in individual channels. 

In this Letter, we consider the scattering of identical 
atoms from a spherically symmetric, weakly-interacting 
atomic Bose-Einstein condensate (BEC), held in a finite, 
localised trapping potential. In this situation, the phase 
shift 6 (k) of /-wave scattering at momentum k = can 
be related to the number n c of bound collective excita- 
tions of the condensate by 



6 l (0)=ir(n c + 6 l0 + a/2), 



(1) 



where a — 1 for a bound state exactly at the threshold 
for s-wave scattering and a — otherwise. Equation 
is the generalisation of Levinson's theorem to the Bo- 
goliubov equations describing the excitations of a weakly 
interacting BEC and is the main result of this Letter. 

Scattering of cold atoms is a fundamental physical pro- 
cess relevant to high-precision atomic spectroscopy and 
quantum information processing. As temperatures are 
lowered, condensation of bosonic atoms is inevitable and 
the scattering of single atoms from condensates needs 
to be understood. Scattering experiments involving con- 
densates have already been demonstrated in the context 
of four- wave mixing experiments 0, 0. With the de- 
velopment of atomic lasers [Ej as coherent matter wave 
sources and with the flexibility introduced by trapping 
and guiding of cold atoms with micro-fabricated electrical 
circuits precision measurement of scattering prop- 

erties becomes feasible. Moreover, interferometric mea- 
surements should allow a direct access of the phase shifts 
and thereby confirm Levinson's theorem experimentally 
in contrast to the case of conventional atomic scattering 
experiments where usually intensities are measured only. 



Theoretical attention has been given to identical particle 
scattering from BECs at low energy where transparency 
effects have been predicted || and at high energy where 
density distributions and quantum correlations can 
be probed [lfj]. Very recently, negative time delays in 
one-dimensional scattering from atomic BECs have been 
predicted 

EH 

In the following, we will apply a multi-channel scat- 
tering formalism to the Bogoliubov equations and prove 
the relation by contour integration. We will discuss 
special situations that can occur like unstable complex- 
energy collective modes and continuum bound states. An 
instructive numerical example of a realistic scattering sit- 
uation is given and the role of the condensate and channel 
coupling for the Levinson theorem are discussed. 

In the standard Bogoliubov approach , the weakly 
interacting Bose gas is described by a condensate with a 
small amount of coherent quantum depletion and a gas of 
non-interacting quasi-particles describing small thermal 
or externally induced collective excitations. The quasi- 
particles are mixtures of particle and hole excitations as 
long as they are located in the condensate region. Out- 
side they just become free particles and can be identified 
with elastically scattered identical atoms. The Bogoli- 
ubov picture relies on the number N of condensate atoms 
being large; indeed, the rate of inelastic scattering on a 
BEC decays exponentially in the Born approximation as 
N grows large |tj ■ 

The quasi-particle energies t v i and the particle (it) and 
hole (v) amplitudes are obtained by solving the Bogoli- 
ubov equations for each partial wave / 

(Tj + Vii)u„j(r) - V 12 v„i(r) = (e vJ + n)u„i(r) (2a) 
(Ti + V 2 2)v vl (r) ~ V 2 iu vl {r) = {-e vl + n)v vl (r), (2b) 

where T t = -Ti 2 /{2m)d rr + 1(1 + l)/(2mr 2 ) is the kinetic 
energy with the atomic mass m. The potential is the 
same in both channels Vn = V22 = VtrapM + 2gn(r) 
and so is the off-diagonal coupling V12 = V21 = gn(r). 
The condensate density n(r) and the chemical potential 
fi are determined by solutions of the stationary radial 
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Gross-Pitaevskii equation: 

{-^drr + W) + ^^r} <P<r) = W (r) , (3) 

where g = 4irh 2 a s /m, a s is the s-wave scattering length, 
and the condensate order parameter ip(r) is assumed 
real and normalised to one. Here we make use of the 
scattering length approximation for the interatomic in- 
teraction and assume zero temperature, where the den- 
sity n(r) = Nip(r) 2 /(Airr 2 ). Both these assumptions are 
taken for simplicity. Generalisation to finite tempera- 
ture Hartree-Fock-Bogoliubov schemes and a more elab- 
orate treatment of interactions are possible and straight- 
forward as long as the structure of Eqs. J2J i s preserved. 
The trapping potential Vt ra , p (r) is a finite localised well 
that falls off to zero sufficiently quickly for large r 0] 
and traps the condensate, leading to a negative chemi- 
cal potential fj, < 0. Formally, the Bogoliubov equations 
© can be understood as scattering equations for two 
coupled channels u and v. At e„z + /i > there is a 
scattering continuum where the u channel is open and 
the v channel is closed. At large distance r, therefore, 
the hole amplitude v u i{r) decays exponentially while the 
asymptotic form 

u u i — ► sin[fcr + Itz/2 + S l (k)} as r — > oo (4) 

of the particle amplitude defines the phase shift 
S (k) of the scattered atom with wavenumber k = 

J 2m/Ti 2 (t v i + fi). Note that there is a second scatter- 
ing continuum for — e„z + fi > where the v channel is 
open. It is related to the first continuum by the general 
symmetry of the Bogoliubov equations J5J that allows 
one to construct new solutions by interchanging u and v 
and simultaneously changing the sign of e. This property 
can be traced to the invariance of the coupling matrix V 
under exchange of diagonally opposite matrix elements 

V12 «-> V21 and V11 <-> V 22 . (5) 

We now proceed with the proof of Eq. 0) for the spe- 
cial case 1 = 0. To this end, we slightly generalise the 
Bogoliubov equations J2J by introducing the channel mo- 
menta k\ and k 2 as independent variables and write in 
matrix notation 

{-d rr t + V)4> = K 2 cj>, (6) 

wherein = diag(Jbi, fo) and = \u v {r), v v (r)] T . We have 
now chosen units where h 2 /(2m) = 1 and dropped the 
channel index I. Analytical properties of multi-channel 
scattering solutions are most conveniently discussed by 
analysing the Fredholm determinant A(fci,fc2), which 
generalises the Jost function of single-channel scattering 
[2(. We define the Fredholm determinant by A(ki, fe) = 
det T with the Jost matrix T 

/•OO 

T = \+\ e iKr V{r)^(K,r)dr. (7) 
Jo 



Here, the columns of the matrix r) are regular so- 
lutions of Eq. (JHJ. Q(K, r) obeys the system of coupled 
integral equations 

K$(K, r) = siTi{Kr)+ [ sm[K{r-r')]V(r')<f>(K, r') dr' . 
Jo 

( 

The Fredholm determinant A(fci, fe) is an entire function 
of k\ and k<z for finite range couplings. Under weaker 
conditions [l9| the function A(fci,fc2) is still analytic for 
3*i > and 3/fc 2 > 0. 

The function A(fei, fe) contains the complete informa- 
tion about the physical scattering process and allows one 
to construct the S matrix. In particular, if k\ or ki are 
real 

A(fe 1 ,fc 2 ) = |A|e- ia (9) 

is directly related to the phase shift 5 of the open channel. 
Furthermore, we have the symmetry 

A(k 1 ,k 2 ) = A*(-k* 1 ,-k 2 *), (10) 

which follows for a real symmetric coupling matrix Vij, 
and for Sfej > and 3*2 > 0, the relations 

A(fci,*a)->1 as |fei|,|*3|-f 00 (11) 

hold 0. The symmetry 

A(k u k 2 ) = A(fc 2 ,fci) (12) 

is a special property of the Bogoliubov equations and will 
be instrumental for the proof of Eq. QJ. Equation |Q 
can be derived from the symmetry JSJ using Eqs. J7J and 

©■ 

We now return to the Bogoliubov equations J5J where 
the channel momenta k\ and ki are not independent vari- 
ables but related through k\ + k\ = 2/i. Resolving for 
hi, we specifically choose 

k 2 {k 1 )=i^k 2 -2 l i, (13) 

where \f ■ denotes the usual positive square root with a 
branch cut along the negative real axis. If k 2 is given 
by (|13fl we say that fci is in the physical sheet. The 
choice (fT3)l assures that 3*2 > 0, whenever 3*i > 0. We 
see that A(k±) = A(fci, k 2 (ki)) is analytic in the upper 
complex half plane, except for a branch point at ki = 
i\J — 2fi and a cut along the imaginary axis where k\ji > 
y/— 2/i. originating from the square root (|13|) as shown in 
Fig-H This is exactly where k 2 becomes positive, i.e. the 
v channel is open and supports scattering solutions. 

Bound solutions of the Bogoliubov equations corre- 
spond to zeros of A(fci) in the upper half of the phys- 
ical sheet where 3*i > 0. Solutions with positive quasi- 
particle energies e v , the usual case for stable ground 
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FIG. 1: Analytic structure of A(fci) in the upper half of 
the physical sheet and integration contour C = S*=i ^ as 
detailed in the text. The segments C7, Ci and C3, C5 lie exactly 
on the real and imaginary axis, respectively. 



state condensates, are found on the imaginary axis at 
ki/i € [0, y/— fi\. Due to the symmetry JSJ of the Bogoli- 
ubov equations, every bound state at k±/i = yj—e u — [i < 
has an image at k\ji — y / +e L , — /i > y/—fi. Phys- 
ically, only the solution with positive quasiparticle en- 
ergy is meaningful. At k\ = i^—fx corresponding to zero 
quasiparticle energy eo = 0, a trivial bound solution of 
the Bogoliubov equations causes a double zero of A(fci). 
This solution with uo = vq oc <p(r) is proportional to the 
condensate order parameter and does not describe a con- 
densate excitation. In fact, quasiparticle excitations are 
confined to the orthogonal complement of <p(r). 

following integral in the complex fci-plane over the 
closed contour C = J2i=i shown in Fig. 

J d\n A(fci) = 2nin z = 27ri(2n c + 2) , (14) 

which counts the number n z of zeros of A(fci) with their 
multiplicity within the region enclosed by C. From the 
preceeding discussion it is clear, that n z is given by twice 
the number of bound collective excitations n c enclosed by 
C plus the contribution from the trivial solution. Note 
that the contour C is more complicated than in the proof 
of the ordinary Levinson theorem, due to the cut and the 
branch point which has to be avoided. 

We will now discuss contributions to the integral on the 
left hand side of Eq. (|14[) . As in the usual proof of the 
Levinson theorem, the contributions on the negative and 
positive imaginary axis can be related to the scattering 
phase shift S and with Eqs. © and (|10|l we obtain 

f dlnA(fci) =2i [5(e) - S(R)] . (15) 

The small semicircle Cg of radius e around the origin gives 
a vanishing contribution for small e as do the quarter 
circles C2 and Cq for large R due to Eq. lfl"TJl . 

The remaining contributions of the integrals over C3, 
C4, and C5 along the imaginary axis can be related to 



the contributions along the real axis due to the special 
structure of the Bogoliubov equations. In fact, Eq. ifLU) 
defines a conformal mapping, which maps the segments 
C3, C4,, and C5 of the ki plane onto the segments C7, Cg, 
and C\ of the &2 plane, respectively. With the symmetry 
i[T2)l we obtain J Ca din A(fci) = J^dlnA^i) and like- 

wise Jc 4 =Sc a and /c 5 =/ Cl - 

We have now evaluated the left hand side of Eq. 114(1 , 
which becomes 4i[<5(0) — 5(oo)] in the limits e — > + and 
R — > 00. Note that we can set 6 (00) = due to Eq. lllfl . 
Thus we obtain Eq. for I = and under the assump- 
tion that there are no zeros of A(fci) on the real axis. 
The generalisation of this proof for I > is straightfor- 
ward within the partial-wave formalism. However, for 
I > all zeros enclosed by the contour C correspond to 
collective excitations since the modes u and v are always 
orthogonal to the condensate wavefunction <p. 

Finally, we consider the possibility of zeros of A(fci) 
on the real axis. Additionally to zeros at the threshold 
k\ = 0, also encountered in single-channel scattering, we 
cannot exclude the possibility of zeros of A(fei) at real 
k\ 7^ 0. If zeros of A(fci) on the real axes exist, they fully 
contribute to a continuously measured phase shift in the 
same way as continuum bound states do in conventional 
multi-channel scattering theory with the exception of a 
zero at k\ = in the I = channel 2] . This specific case 
is known as a half bound state as it only contributes 7r/2 
to the phase shift, hence the term a/2 in Eq. 0J. 

So far we have assumed that the condensate is in a sta- 
ble ground state of the trap, in which case all bound exci- 
tations described by the Bogoliubov equations have finite 
pseudo norm 77 = J (\u\ 2 — |u| 2 ) dr ^ and those with 
77 > have positive quasi-particle energy e v . However, 
this restriction is not necessary. For excited states ip{r) of 
the stationary Gross-Pitaevskii equation ©, anomalous 
solutions of the Bogoliubov equation are possible with 
e„ < and 77 > 0. Such anomalous modes are known 
to occur for solitons in highly elongated traps 0, IHj 
and for vortices [l5j . Furthermore, bound solutions with 
complex e v and r\ = may occur, which correspond to 
zeros of A(fci) on the physical sheet but off the imagi- 
nary axis. These solutions arise when the condensate is 
in a stationary but unstable excited state of the Gross- 
Pitaevskii equation and describe modes of exponential 
decay, predicted and seen, e.g. , for dark solitons [lH Ht| 
and attractive condensates |l8| . These unstable complex 
modes are not to be confused with scattering resonances, 
which have zeros below the real axis and in the unphysical 
sheet. Both anomalous and complex Bogoliubov modes 
contribute to the contour integral (|14f) and thus are pre- 
dicted to be visible in the build up of the elastic scatter- 
ing phase shift (5(0), which can, in principle, be probed 
experimentally. 

This is exemplified in Fig. [21 for a ground-state BEC of 
2000 Rubidium atoms. Note that the relevant range of 
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FIG. 2: Phase shifts Si from the full solution of Eqs. © 
(solid), simple potential scattering without coupling [V12 = 
in Eqs. @] (dashed). In the example, 2000 Rubidium atoms 
are held in a finite harmonic trap of frequency lo — 27r200s _1 
and depth Vb = S.Ohui resulting in the chemical potential 
H — — 3.34hui and allowing for up to two bound collective 
excitations for 1 < Z < 5. 



is based on a weakly interacting Bose gas without inelas- 
tic scattering processes. It should be seen as a challenge 
both to experiments and theory to find the corrections 
to Eq. Q in a real interacting system. Dilute-gas atomic 
BECs are an ideal system for this because the experi- 
mental set up is stupendously manageable. 




r [v /3] 



FIG. 3: Number of bound collective modes for a square-well 
of depth Vb = 10, 100, 1000 h 2 /ma 2 , respectively, and radius 
a with a Thomas-Fermi approximation for the wave function 
as a function of the nonlinear coupling V = Na 3 /a. Full 
solution (symbols), potential scattering (dashed) as in Fig. [5] 



scattering energies corresponding to temperatures of the 
order of 10 nK is within the accessible range of current 
cold-atom experiments. For I — one bound collective 
excitation and the trivial solution account for a phase 
shift at threshold <5 (0) of 2ir both in the full Bogoliubov 
equations (solid) and in the potential scattering approx- 
imation (dashed) where we have set V12 = 0. Both solu- 
tions coincide in the absence of a bound mode for 5i(ki) 
which exhibits, instead, a resonance caused by a zero of 
A(fci) with fci having a small negative imaginary part. 
The interesting case is #2(^1) which indicates one more 
bound state in the full solution compared to simple po- 
tential scattering. 

The offdiagonal coupling V12 is of attractive nature and 
binds extra states although the interparticle interaction is 
repulsive in both examples. For large scattering momenta 
and large angular momenta this coupling to the v mode 
becomes less important, as the details of the interaction 
region are hardly probed. Figure [3] shows the number of 
bound states of a square well model system. As the well is 
filled up with repulsive atoms, fewer collective excitations 
are bound. Note that for T > Vb/6 the trapped BEC 
has bound quasiparticle states only due to the coupling 
between hole and particle modes. 

We have presented a Levinson theorem relating the 
number of bound collectively excited states of an inter- 
acting many-body system, a BEC, to the phase shifts of 
single-particle scattering. With the possibilities of cold 
atom scattering and interference, we can expect to see 
a direct experimental verification of a fundamental the- 
orem of mathematical physics. However, our derivation 
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